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Abstract 

We obtain a new solution to the star-triangle relation for an Ising-type model with two 
kinds of spin variables at each lattice site, taking continuous real values and arbitrary integer 
values, respectively. The Boltzmann weights are manifestly real and positive. They are expressed 
through the Euler gamma function and depend on sums and differences of spins at the ends of 
the edge. 
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The star-triangle relation is a distinguished form of the Yang-Baxter equation for Ising-type 
models on two-dimensional lattices. In these models the fluctuating variables or "spins" are as- 
signed to lattice sites, while two spins interact only if they are connected by an edge of the lattice. 
Remarkably, many physically interesting models in this class can be solved exactly, for instance, 
the 2-d Ising [Ij and the chiral Potts [2l[3] models (see [UlS] for a review of other known cases). The 
star-triangle relation plays the role of the integrability condition for these models. 

Recently, Bazhanov and Sergeev found a very important "master" solution [6] of the star- 
triangle relation, which contains all previously known solutions of this relation as particular cases, 
and provides many new interesting examples. The above master solution is expressed in terms 
of the elliptic gamma function. It contains two arbitrary free parameters p and q, which play 
the role of elliptic nomes. The spin variables take continuous real values on the circle. Various 
interesting cases arise when the parameters p and q approach some special values. These limits 
are rather singular and their consideration requires a detailed analysis on a case by case basis (for 
a comprehensive account see the forthcoming paper [5]). In the present paper we consider one 
particular remarkable limit of the master solution, which leads to a model with two types of spin 
variables at each lattice site. These variables take continuous real values and discrete integer values, 
respectively. Mathematically the master solution of |6] is related to the elliptic beta integral [7]. Here 
we consider first a degeneration [8] of this solution related to the "hyperbolic beta integral" [9|[T0]. 
Then we consider the so-called "strong coupling regime", where new integer- valued spin variables 
dynamically arise at each lattice site. Our calculations generalise the results of 

Consider the square lattice of N sites and assign spin variables Xj, j = 1,2, . . . , N , taking some 
continuous set of values, to all lattice sites. Two spins interact only if they are connected by an 
edge. Let Wa{x,y), yVa{x,y) denote Boltzmann weights for horizontal and vertical edges, where 
X and y are spins at the end of the edge, as shown in Figure [TJ We assume that the edge weights 
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Figure 1: Edges of the first (left) and second types and their Boltzmann weights. 



depend on the (additive) spectral variable a and are related to each other by the crossing symmetry 
yVaix,y) = yVrj-aix,y), where r] is the model dependent crossing parameter. Moreover we assume 
that the weights are reflection symmetric Wa{x,y) = Waiy,x). It is also convenient to introduce 
single-spin weights S{xj), which are independent of the spectral variable a. The partition function 
is defined as 

" n^°(^*'^j) n W,,_„(a;fe,xj) Jj5(xm)iixm, (1) 

(ij) (kl) 

where the first product is taken over all horizontal edges (ij), the second over all vertical edges 
(kl) and the third over all internal sites of the lattice. We will implicitly assume fixed boundary 
conditions. Obviously the single spin weights S{xj) can be included into the definition of the edge 
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weights, but we prefer not to do so. The model is integrable if the weights satisfy the star-triangle 
equation of the form 

dxoS{xo)Wrj-ai (xi, Xo)Wr,-a2 (a;2, Xo)>V,,_q3 (xs, Xq) 

where the spectral parameters ai, 02, as satisfy the relation 01 + 02 + 0:3 = 77 and the factor 
7^(ai, 02, as) is independent of the spins xi,X2,xs. 

Let us now define the Bazhanov-Sergeev master solution [6j. Introduce the elliptic nomes p,q, 
and the crossing parameter rj 

ITT 

q = e--, p = e-% |p|,|q|<l, e'^'? = pq , r? = --(r + a). (3) 
The elliptic gamma function [12l[13] is defined as 



-2\nx 



°° 1 _ e2i2q2i+ip2fc+i g- 

^^^^ " l-e-2iV^+lp2fc+l " """"P I „(pn_p-n)(qn_q-n) \ ' 

where the exponential form is valid in the strip | Imz \ < Rer/. The Boltzmann weights are defined 
by 

^/ N e''/^ ,,,, , , . ,^(x + y + ia)<l>(x - y + \o) 

S{x) = —M2x p)M^x q) , Wa{x,y) = k{o)- \) J . ^— ^ , 5 

ZTT + y — iaj<P(x — y — \o) 

where 'di{z \ q) is the standard Jacobi theta function [14J. Note, that the weights are periodic 

W„(x + 7r,y) = >V„(x,y + 7r) = >V„(x,y), 5(x + vr) = (6) 

Correspondingly, the spin variables in ([1]) take values in the interval < Xm < tt. The normalisation 
factor K{a) in ([5|) reads 

log K(a) = — . (7) 

^n(p"-p-")(q"-q-")(p'"q" + p-"q-") 

The above Boltzmann weights satisfy the star-triangle relation ([2|) where the integration is taken 
over the segment < xq < vr. Below we assume that the spectral parameter a lies in the domain 
< a < r/. In this case the weights ^ are real and positive. 

With the normalisation ([7]) the factor 7^ in ([2]) is equal to one. For the same normalisation, the 
weights ([5]) also satisfy the following inversion relations 

Wa{x,y)W-a{x,y) = 1, 

1 (8) 

dzS{z)Wri-a{x,z)W.,^+aiz,y) = {6{x + y) + 6{x - y)) . 

2S[x) 

These relations allow one to show [6] that in the thermodynamic limit, when the number of lattice 
sites goes to infinity N oo, the bulk free energy of the system vanishes, 

lim iV"^logZ = 0. (9) 

A^— >-oo 
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As a mathematical identity the star-triangle relation ^ is equivalent to the elliptic beta integral 
[7j, which plays a fundamental role in the theory of elliptic hypergeometric functions. Its significance 
for solvable models of statistical mechanics was discovered in [6] . 

Now consider the hyperbolic limit [8j of Introduce the modular parameter b to take values 
in one of the following two regimes 

(i)b>0, (ii) \b\ = l, lm(b^)>0. (10) 

The non-compact quantum dilogarithm is defined as the integral 

(^(.)=exp{i/ . J'T-"^^ ( /h J ' |Im(^)| < Re(ry), (11) 
I 4 Jr+\o y smh(yb) smh(y/b) J 

where the singularity at the origin is taken below the contour. Substituting 

p = e"^^ , q = e"'' , Xi = x^e , = a^e , i = 1, 2, 3 , (12) 
into ([5]) and taking the limit e — )• 0, one obtains [8], 

S{x) = 2 sinh(27ra;ob) sinh(27rxo/b) , >V„(x, y) = /t(a)-^e^""'^ + y + la) ip{x - y + la) ^ 

ip{x + y — \a) ip{x — y — \a) 

These weights satisfy the star-triangle relation ([2]), where the integration taken over the whole real 
line and the crossing parameter is given hy rj = {h + b~^)/2. The normalisation factor K(a) in (jl3p 
reads 

If e'^y^dy , 2 T^'^ 2 ^'^ 

logK(a) = — / T-n r-, — jn r -|- ttiq r? H . (14) 

8yR+ioysinh(yb)sinh(y/b)cosh(2OT) ^ 3 ' ^ 24 ^ ' 



The weights ()13p define a physical model of statistical mechanics, rather similar to the Faddeev- 
Volkov model [15tll6|. The spins now take arbitrary values on the real line Xj E M. The above 
limiting procedure does not affect the bulk free energy of the system, so that the relation ([9]) 
continues to hold for the new model. 

Let us now consider a further specialisation of the above generalised Faddeev-Volkov model ()13p . 
when the parameter b — )■ i, or equivalently, ry ^ 0. In this limit the weights Wa(x,y) develop a 
series of delta function like peaks for integer values of the sum x -|- y € Z or differences x — y € Z of 
the arguments x,y. This limit was previously investigated for the Faddeev-Volkov model [TT] and a 
similar phenomenon was observed. A correct limiting procedure capturing a fine structure of these 
peaks requires a redefinition of the spin variables. The asymptotics of the non-compact quantum 
dilogarithm ip{z) and the function K{a) read 

^{n + xr?) = e--/i2(4^r?) - , = {87rrj) , (15) 

where r/ — )• 0, n G Z, x <C rj^^ and < /3 < 1. Using this in ()13p one obtains the following — )• 
asymptoticsjj 

S{n + XT]) -^87r^(x^ + n^) rj^ , 

2-5/3 r(^~^ ± _[_ i(a:-y)-(m-n) -^ (Ig) 



y)+(m— n) n 



^Here we use the compact notation, where the ± symbol in the argument of the gamma function indicates the 
product of gamma functions with both signs are taken as e.g. V(a ±b) = T{a + b)T[a — b) 
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where m,n £ Z and x, y € M. Substituting this into ^ and cancehing out singular factors from 
both sides, one obtains the following star-triangle relation 



noeZ 



(17) 



= W/3i(x2,n2 I X3,n3)WfS2{xi,ni \ X3, n3)W/33(a;2, n2 | xi,ni) , 
where /?i + /32 + /^a = 1, and the Boltzmann weights are given by 

= — (2;2 + n2) , n | y, - ^ 2 ^ ^ 2 2 ^ I 2 2 



27r PVi ; r(lz^) r(^+'^ ifc '(^+?/)+(m+n) ^^p^l+/3 _^ \{x-y)+{m-n) ^ 

(18) 

These weights give a new solution of the star-triangle relation for an Ising type model, where sites 
of the lattice are assigned "dual spin" variables x 

X = {x,n) , X € M , n G Z , (19) 



taking arbitrary real and integer values, respectively. The Boltzmann weights (jlSp satisfy the spin 
reflection symmetry x -H- y, for two spins x, y connected by an edge as in Figure [TJ The bulk free 
energy (j9|) of the corresponding 2-d lattice model remains unchanged. 

The star-triangle relation ()17p is derived in the strong-coupling limit of the generalised Faddeev- 
Volkov model (|13p . the latter model arising as a hyperbolic degeneration of the master solution 
([5]). From either of the models ([5]), (I13p . new solutions of the star-triangle relation can be found 
as different singular limits of the parameters p,q, and b respectively. In each case, the spins and 
spectral parameters should be redefined, as was done here in (fT2]) and (fT6]l . in order to obtain 
convergent limits. 
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